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Abstract 

We consider an American put option under the CEV process. This corre- 
sponds to a free boundary problem for a PDE. We show that this free boundary 
satisfies a nonlinear integral equation, and analyze it in the limit of small p = 
2r/a 2 , where r is the interest rate and a is the volatility. We use perturbation 
methods to find that the free boundary behaves differently for five ranges of time 
to expiry. 



1 Introduction 



The pricing and hedging of options has its origins in the Nobel prize winning work 
of Black, Scholes, and Merton [2], who assume that the price of an underlying asset 
S(t) follows a geometric Brownian motion with constant volatility. The price C(S,t) 
of a European call option at time t for an asset with price S, strike K, and expiry 
T is then readily established, and is presented in terms of the normal distribution 
function. However, there is sufficient empirical evidence [I] to suggest that in many 
cases the assumption of constant volatility does not match well to the observed market 
data. Rather, evidence points out that the implied volatility, which is obtained by 
equating the model price of an option to its market price and solving for the unknown 
volatility parameter, varies with the strike price across a wide range of markets. This 
phenomenon is known as the volatility smile or frown, depending on the shape of the 
curve, and is not captured by the Blacks-Scholes model with a constant volatility. 
As a result, there have been various ideas as to how to modify and extend the basic 
Black-Scholes framework, to account for this phenomenon. One of these is the constant 
elasticity of variance (CEV) diffusion model, which was introduced by Cox and Ross |3J 
in the context of European options. Unlike Black-Scholes, the CEV model is capable 
of reproducing the volatility smile. 

Other work on European options under a CEV process include Davydov and Linet- 
sky |3], Hu and Knessl [8] and Lo, et. al. |12j . However, there exists little or no 
analytic work for the valuation of American options under a CEV process. The analy- 
sis of these options are more difficult than the corresponding European options in that 
the American options may be exercised prior to the expiration dates. Mathematically 
the American options lead to partial differential equations (PDE) with free boundaries, 
which can only rarely be solved exactly. In this paper, we apply asymptotic analysis 
to a CEV model to examine the behavior of the free boundary under different scaling 
regimes for the time to expiry, in the limit of small p = 2r/a 2 , where r is the interest 
rate, and a is the volatility. This limit has a small interest rate and/or large volatility, 
and is of particular relevance to the financial status of the current economy. We will 
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employ singular perturbation methods, including matched asymptotic expansions. The 
main result is the derivation of a nonlinear integral equation that is satisfied by the 
free boundary, from which we shall analyze its asymptotic structure for five different 
ranges of time. The main results are summarized in section 2 and derived in section 
section 3. 

Asymptotic analysis and singular perturbation methods have been recently em- 
ployed in the context of both European and American options, and this work includes 
Knessl [§1 HO], Howison [7], Kuske and Keller [UJ, Addison, et al. jT], Evans, et al. [5], 
Fouque, et al. [6J, and Widdicks, et al. [13] . 



2 Problem Statement and Summary of Results 

We let P(S, Tq) denote the price of an American put option for an asset with price S at 
some time T prior to expiry Tp. We assume that S satisfies the stochastic differential 
equation 

dS = nS dt + aVS dW t . (1) 

where W t is a standard Brownian motion, a is the volatility of the underlying asset, 
and p = r is the risk-free interest rate. We note that unlike Black-Scholes, this model 
only guarantees non-negativity of S (S > 0), so the chance of absorption at 0, i.e., 
bankruptcy, occurs with positive probability. 
Introducing the new variables 

cr 2 2r 
t=—(T F -T ), P =- 2 , (2) 

we find that P satisfies the following boundary value problem 

P t = SP SS + pSPs - pP; t > 0, S > a(t) (3) 

P(S, 0) = max(iT - S, 0) (4) 

P(a(t),t) = K -a(t) (5) 

P s (a(t),t) = -1 (6) 
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P(0,t) = (7) 

where a(t) is the free boundary in the new time variable. We also have P(S,t) = K — S 
for < S < a(t), and a(0) = K. For S < a(t) the option should be exercised, and for 
S > a(t) it should be held. 

We convert ^ — ^ into an integral equation by first making a change in coordi- 
nates, letting 

P(S, t) = K-S + P(V, t), V = S - a(t) (8) 
where V > 0. Then P satisfies the PDE 

P t - a'(t)P v = [V + a(t)}P vv + p[V + a(t)]P v - pK - pP; V, t > (9) 



with the initial and boundary conditions 

P(V,0) = V 

P(0,t) = Fv(0,t) = 0. 
We introduce the Laplace transform 

POO 

Q(9,t) = / e- 0V P(V, t) dV. 
Jo 

Using (12) in ^ and (11) then yields 



(10) 



'111 



(12) 



Qt + (e 2 + P 6)Q e = [a'(t)6 + (6 2 + p6)a(t) - (26 + 2p)]Q - 



pK 



with the initial condition 



Q(6,0) 



e 2 ' 



(13) 



(14) 



Using the method of characteristics, it can be shown that the only acceptable solution 

'0 + P z 



to (13) is 



Q(9,t) 



(t)e 



e 2 



e/ P 



z + 1 



cxp 



pza ( t + p 1 lo£ 



e z + i 



dz. (15) 



The next result readily follows. 
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Theorem 1. The option price P(S,t) for the CEV model has the integral representa- 
tion 



P(S,t) = K 



e ev Q(9,t) d9, 

2m J Br 



(16) 



where 31(9) > on the Bromwich contour, and Q(9,t) is given by (15). 



Moreover, after setting t = and using a(0) = K and (14) in (15), it follows that 
a(t; p) satisfies the nonlinear integral equation (IE): 



,-A"0 



Kp 



z + 1 



cxp 



-pza (^p 1 log ^ 



9 + p z 
9 z + 1 



dz. 



(17) 



In the next section we use asymptotic methods to analyze this IE for five different 
scales of time t, in the limit of small p. We let p = e~ x so that A = — logp — > oo. The 
final results for the free boundary a(t; p) are listed below, and we sketch the derivations 
in section 3. 



(i) t = u/\ = 0(A _1 ), < u < K: 

/x , i— rr7x 9 log Xtu — y/ Koo 1 V 'Kuj — u , 
a(t;p) = (v^-VK) 2 +^ - + - log 



k- Vk^J 



(18) 



ii) t = K/X + 0(\- 2 ), XH -XK = X(cu - K) = A: 

a(t;p) = ^-F(A), 



where J-"(-) satisfies the nonlinear IE 

-Kv r°° y 

exp 



k Jo e 

For A ±oo, we have 

A 2 A 



-£F -vK 



A 



H*) ~ 4^ + j !og(" A ) " ^ log(87rK 3 ), A ^ -oo 



J*(A) ~ Ae 7 exp 
where 7 is the Euler constant. 



1/ V 



A -)> +00 



(19) 

(20) 

(21) 
(22) 
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(iii) t = u/X = 0(X- r ), K < u < oo: 



(iv) t = 0(1), < t < oo: 

a(t] p) ~ te~ 7 exp 

(v) t = v/p = 0{p- 1 ), v>0: 

a(t; p) ~ -e~ 7 exp 



2 + ^" 



'1 — e u ) exp 



1 



(23) 



(24) 



(25) 



p ~ \e v - I J y ' "~ r \ pK / 

We note that in four of the five cases the expression for a(t; p) is completely explicit, 
and only in case (ii) must we solve a nonlinear IE, which is somewhat simpler than 



the one in (17). We can easily compute P(S,t) as t — > oo, which corresponds to 
the perpetual American option, where the problem reduces to solving an ordinary 
differential equation. Setting P(S, oo) = P°°(S) and using a(oo) to denote the limiting 
value of the free boundary, we obtain from (J3j) — ([^J) 



f°° 1 

P°°(S) = Ke pa{oo) / — e~ zps dz, 
Ji z 2 



where a(oo) satisfies 



Kp / - e - pa ^ z dz = e - pa{oo) . 



For p — > we have 
which is exponentially small. 



1 (-— Mi + o( P )], 



al oo ) = — e ' exp 

p \ pK 



(26) 
(27) 
(28) 



3 Analysis 



3.1 Analysis for t = uj/X, < lu < K 



We first examine (17) on the t = 0(X x ) scale, for small p. Recalling that A = — logp, 
we let 

6 = A/3, z = X ^ + x \ a f t . p ) „ aQ ( w; p) (29) 
P 



for uj 



log p)t = 0(1). Then (17) can be approximated by 



e 

K 



3 A*(x;/3 l(0 ) 



[1 + 0(e~ A )] dx, 



(30) 



/o (3 + x 
where $(x; /5, p)=K/3-(/3 + 

x ) a o{ p(x+p) ' ^ or l ar S e ^ an d fixed /3, we evaluate the 



right hand side of (30) by an implicit form of the Laplace method, assuming for now 



that there is a saddle point where 



3$ 

dx 



-oc 



1 



On 



0(x + p)J x + /3 u \p(x + p) 



X 



0. 



(31) 



Let us denote x = x*(/3) as the solution to (31). It follows that at x = x*, $ ~ 1 so 
that 

l=Kp-(p + x*)a 



P(x* + P) J ' 



Now let u 



f5(x t +/3)- 



Then from (31) we have P 



u>a' (u)— ao(ui) 



(32) 

which we use to 



eliminate P in (32) to obtain the ODE 



Rewriting this as 



[1 - a' Q (u)\[ua' Q (u) - ao(w)] = Ka' (u) 
ua (u) - a Q (u) - 



(33) 



(34) 



1 -«'oH' 

we recognize this as the Clairaut equation. The solutions consist of a one-parameter 
family of lines and the singular solution 



(35) 



which is the envelope of this family. The linear solutions oiq(oj) = ujC — must 
be rejected, since these lead to «o(0) ^ K. The above analysis applies only for < 



u) < K, since the solution (35) vanishes as uj approaches K. Hence we expect different 
asymptotics for uj ~ K. 



We next analyze some higher order terms in the expansion of «o- We evaluate (30) 
by using the Laplace method, which gives 



e 

K 



1 



2vr 



ft + U -\§ xx (x*; p,p) 



[1 + 0(X~ 1 )], 



(36) 



and expand a as 



a (u; p) = a (u) + a^u) + ja 2 (u) + o(A 1 ). 



(37) 



In order to balance the two sides of (36), we need ot\ to cancel the yl/A factor. Hence, 
1 



exp 



-((3 + a;*) (log A) &\ 



/3(z* + (3) 



1. 



Writing (38) in terms of w we obtain 



ol\{lo) = -(cu — V Kcu). 



To find the third order term a 2} we balance the 0(1) terms in (36), so that 

1 



1 



K + 



2vr 



exp 



-(/? + x*)a 2 



/%*+/?) 



It can be shown that $^(2;*; /3, p) ~ — 2 (x* + /3) 2 and then 



«2(w) 



log 



47rir 2 o; \ 



K 



With (35), (37), (39), and (41) we have established (18) 



(3? 



(39) 



(40) 



(41) 



3.2 Analysis for t = uj/X, lu ^ K 



We return to (17) and introduce the scaling 



with 



9 = \p, P = ^ + j, u; = \t = K+j 



a(t;p) = ^(A;p) = - 2 F I A 2 ( t - - ) ; p ) . 



(42) 



(43) 



Then we have e eK = pe Kv '. Also by setting z = (9 + y)/p in (17) this equation 
becomes 



1 r°° 1 

K 



# + 2/ + P 



exp 



dy. (44) 
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With the scaling in (42) and (43) and the fact that p = e is exponentially small, we 



obtain 



X 2 a 



9 9 + y 



+ 0(e 



X 2 a 
T 



K 
T 



1 

A 2 



uK 2 + -\ +o(\- 2 );p 



-vK 2 - 



(45) 



where J-"(A) is the leading term in an expansion of J-"(A; p). Then scaling y = A 2 £ in 



(44) and letting p — > (A — > oo) we obtain the limiting IE in (20). It does not seem 



possible to solve (20) explicitly for .F(A). But we can infer the behavior as A — > — oo by 



evaluating the integral in (20) by an implicit Laplace type expansion, similarly to what 



we did in section 3.1. This will verify the asymptotic matching between the w-scale 



(for oj < K) and the A-scale, and lead to (21). Now consider the limit A — > +oo. For 



\u\ < 0, we rewrite (20) as 



K\v\ 



\u\K 2 V 



cxp 



--mu\K 2 - V ) 



dr) 



+ / - <; exp 
/ u 

roc g-t 

\v\K 2 



I T(\u\K 2 (l-u)) 
u 



W\K 2 



cxp 





\K 2 ) 


u 


v\ 


K 2 



(46) 



dv. 



Here we broke up the integral over (0, oo) into the two ranges (0, |^|i^ 2 ) and (|z/|i^ 2 , oo) 
and made some elementary substitutions. Now, for A — > — oo we have .F(A) ~ so 



that the first integral in the right hand side of (46) will vanish as v — > — oo. If 



J-"(A) — > as A — > +oo the second integral in (46) will also vanish, and the third may 
be approximated by using 



oo e -v 



dv 



loge - 7 + 0(e), e-^0 + . 



Hence (46) can be replaced by the asymptotic relation 

p K\v\ 

-log 



K 





\K 2 ) 




p\ 


K 2 



-1 



(47) 



(48) 



which upon exponentiation leads to the asymptotic result given in (22), for ^(A) as 
A -)■ oo. 
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3.3 Analysis for t = cu/X, K < u < oo 

In the remaining time ranges, a(t; p) will be exponentially small as p = e~ x — > 0, and 



our analysis of (17) will rely heavily on the asymptotic form in (47). We let z = Z/p 

e + P z 



in (17) to obtain 

„-K6 



Kp 



Z + p 



cxp 



Za I p 1 lo 



9 Z + p 



dZ. 



(49) 



Now we scale Z = Az* and 9 = X9 st , let a(t; p) = a(Xt; p) and note that in sections 3.1 
and 3.2 we have already characterized a(Xt; p) for Xt = u < K and u ~ K. We also 



simplify the argument of a(-) in (49) using 

P 



a 



1 

P 1 



a 



1 1 

9 ~ Z 



0{ P ) 



a 



1 

Z. M 



0(e- x X) 
(50) 



When Z = 9 we have z* = 9* and we rewrite the integral in (49) by splitting the range 
of integration into z* G #*/[! — K9*\) and z* G (#*/[! — K9*\, oo), thus obtaining 



-K\e, 



Kp 



1-K6t, 



+ 



9* 



2# 



exp 



-Xz*a 



1 



fife*. 



(51) 



In the first range a{u) ~ (yK — ^/u) 2 and the first integral will be o(l) as A — )■ oo, since 
9~ l — z^ 1 < K when z* < 9*/ [I — K9*\. In the second integral a will be exponentially 
small and the main contribution will come from very large values of 2*, where roughly 
= 0(<5: -1 ). Then we write a{9~ 1 — z^ 1 ) ~ o-iO' 1 ) and using (47) we conclude that 



Kp 



1 



1-K8* 



cxp 



-Xz*a 



1 



dz* ~ — loe 



Aa 



1 



- 7 - lo£ 



1-K6* 



(52) 



with an error that is o(l) as A — > oo. Then exponentiating (52) and replacing 9* by 



XK 



K- 



we obtain the asymptotic result in (23). 
For to — y K we note that p w l K ~ x = p~ 1 e~ K / t = p _1 exp 
= p _1 exp [—A + ^ + 0(A -1 )] ~ exp (A) and (u — K)/X = A/A 2 , which can be used 
to verify the asymptotic matching between the A-scale and the w-scale for u > K, in 
the intermediate limit where u> j. K and A — > oo. 
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3.4 Analysis for t = 0(1), < t < oo 

Next we consider times t = 0(1). We scale z = 8w/p. Since we again expect a(t; p) to 
be very small we assume a "WKB-type" ansatz of the form 



"(*; P) ~ 9(t) exp 



(53) 



Expanding a(t; p) in (17) for fixed and p — >■ 0, and noting that 
1 



we have 



—9wa 



log 1 



log 1 + 



log 1 + 



0w 



1 1 ^ ^nf 2 



2 P 



2 / ' 



log 1 + 



9w 



-9wg \ - ] exp 



p V# ^ 2 # 2 



I ) (>X P 



20 2 ' / V0 ' 



Here we also used f(9 1 — (9w) 1 ) ~ f(9 1 ), since w will be scaled to be exponentially 
large. Then setting e = 9g (|) exp — -/ (|) exp [^f (|)], scaling u> = and 



using (47), (17) asymptotically becomes 

-K6 

- 7 - log 



V 1 



(54) 



From the 0(p x ) terms in (54) we conclude that f (1/9) = K l e K and then the 0(1) 
terms determine g(-) from 



e 7 exp 



2^ 2 7 



The above along with (53) establishes the asymptotic result in (24). The asymptotic 
matching between (23) and (24) is immediate, since p K / u '~ l = p _1 e _if// ', and (co — 



K)/ \ ~ uj/ \ = t as uj — > oo. 



3.5 Analysis for t = v/p = 0(p 1 ), v > 

We assume that time to expiry for the option is large, with t = v/p = 0(p^ 1 ). On this 
time scale we assume that 



a(t;p) ~ - exp ( 

p \ pK 



A(v) 



(55) 
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where A(-) will be determined from (17). After scaling 6 = pW, (17) becomes 

P -KpW 



Kp 



w 



z + 1 



exp 



— ze p R A lo. 



W + l z 
W z + 1 



dz. 



(56) 



The major contribution to the integral in (56) will once more come from large values 
of z, so we approximate 



A lo 



W + l z 
W z + 1 



A log 



W 



w 



(57) 



and then applying this to (17) along with (47) leads to 

p -KpW 



Kp 



log < e p k A\ lo 



W + l 

w 



7 



1 

~ ~Kp 

Then after expanding e~ Kpw 



log {(W + l) A (log 



- 7 - \og(W + 1) 
W + l 



w 



+ 0(1) 



(58) 
(59) 



KpW + 0(p 2 ) we conclude that 



A lo 



W + l 
W 



1 



W + l 



(60) 



which determines the function A(-) and establishes (25). 



Finally we verify the asymptotic matching between (24) and (25). For v — > we 
have (e v — 1) _1 = v " 1 — \ + 0{y) and 1 — e~ v ~ v — pt. For t — )■ oo we have 



-K/t 



1 - f + 0(r 2 ) = 1-^p + 0{p 2 ) so that -{\ + \)e~ K l l ~ -v~ l + \ and 



the matching follows. As v — > oo we have A(v) — > e 7 and thus the expansion in (55) 



agrees with the small p expansion of a(oo; p), as given in (28). 
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